home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / plotting / contour / contour.lha / Contour / IOread.c < prev    next >
Encoding:
C/C++ Source or Header  |  1990-08-28  |  11.8 KB  |  375 lines
  1. /*** IOread.c - read input etc ***/
  2.  
  3. #include <stdio.h>
  4. #include <strings.h>
  5. #include "common.h"
  6. #include "contour.h"
  7. #include "plot.h"
  8. #include "Xdefs.h"
  9.  
  10. /* Read the command line */
  11. fileopen(argc,argv) 
  12. int argc;
  13. char **argv;
  14. { 
  15.    FILE   *fp,*optionfile, *fopen();
  16.    double atof();
  17.    extern char   printer[MAXCHAR];
  18.    extern int    cont_out, expand, grid, postscript, oldformat, contlabel;
  19.    extern int    joincurve, printplot;
  20.    extern char   *joinlevel;
  21.    extern double cont_level;
  22.    int    i,FILEFOUND;
  23.  
  24.    /* initialize general X arguments */
  25.    progname     = argv[0];
  26.    font_name    = NULL;
  27.    geometry     = NULL;
  28.    display_name = NULL;
  29.    border_color = NULL;
  30.    fore_color   = NULL;
  31.    back_color   = NULL;
  32.    mouse_color  = NULL;
  33.    icon_color   = NULL;
  34.    reverse      = 0;
  35.    border_width = 0;
  36.  
  37.    FILEFOUND = FALSE;
  38.    for (i=1; i < argc; i++) {                     /* Parse Line */
  39.       if (argv[i][0] == '=') {
  40.          geometry = argv[i];
  41.          continue;
  42.          }
  43.       if (strncmp(argv[i],"-geom",5)==0) {        /* new geometry */
  44.          if (++i >= argc) usage(progname);
  45.          geometry = argv[i];
  46.          continue;
  47.          }
  48.       if (index(argv[i],':')!=NULL) {             /* host:display */
  49.          strncpy(display_name,argv[i],sizeof(display_name));
  50.          continue;
  51.          }
  52.       if (!strncmp(argv[i],"-d",2)) {             /* new display */
  53.          if (++i >= argc) usage(progname);
  54.          display_name = argv[i];
  55.          continue;
  56.          }
  57.       if (strcmp(argv[i],"-rv") == 0) {           /* reverse display */
  58.          reverse = 1;
  59.          continue;
  60.          }
  61.       if (strcmp(argv[i],"-bw") == 0) {           /* border width */
  62.          if (++i >= argc) usage(progname);
  63.          border_width = atoi(argv[i]);
  64.          continue;
  65.          }
  66.       if (strcmp(argv[i],"-bd") == 0) {           /* border color */
  67.          if (++i >= argc) usage(progname);
  68.          border_color = argv[i];
  69.          continue;
  70.          }
  71.       if (strcmp(argv[i],"-fg") == 0) {           /* foreground color */
  72.          if (++i >= argc) usage(progname);
  73.          fore_color = argv[i];
  74.          continue;
  75.          }
  76.       if (strcmp(argv[i],"-bg") == 0) {           /* background color */
  77.          if (++i >= argc) usage(progname);
  78.          back_color = argv[i];
  79.          continue;
  80.          }
  81.       if (strcmp(argv[i],"-fn") == 0) {           /* font type */
  82.          if (++i >= argc) usage(progname);
  83.          font_name  = argv[i];
  84.          continue;
  85.          }
  86.       if (strcmp(argv[i],"-c") == 0) {            /* output contour */  
  87.          if (++i >= argc) usage(progname);
  88.          cont_level = atof(argv[i]);
  89.          cont_out   = ON;
  90.          continue;
  91.          }
  92.       if (strcmp(argv[i],"-e") == 0) {            /* expand into 3D */
  93.          expand = ON;
  94.          fprintf(stdout,"   Commands: expand (3D)     = ON\n");
  95.          continue;
  96.          }
  97.       if (strcmp(argv[i],"-j") == 0) {            /* join curves */
  98.          if (++i >= argc) usage(progname);
  99.          joincurve = ON;
  100.          joinlevel = argv[i];
  101.          if ((strcmp(joinlevel,"HIGH")==0)||(strcmp(joinlevel,"high")==0)||
  102.              (strcmp(joinlevel,"LOW" )==0)||(strcmp(joinlevel,"low" )==0)) {
  103.             fprintf(stdout,"   Commands: joincurve       = ON\n");
  104.             fprintf(stdout,"   Commands: joinlevel       = %s\n",joinlevel);
  105.          } else {
  106.             fprintf(stdout,"   Warning : Invalid keyword -");
  107.             fprintf(stdout," [-j] option only accepts [HIGH/LOW]\n");
  108.             joinlevel = OFF;
  109.          }
  110.          continue;
  111.          }
  112.       if (strcmp(argv[i],"-g") == 0) {            /* grid */
  113.          grid = ON;
  114.          fprintf(stdout,"   Commands: grid            = ON\n");
  115.          continue;
  116.          }
  117.       if (strcmp(argv[i],"-l") == 0) {            /* label */
  118.          contlabel = OFF;
  119.          fprintf(stdout,"   Commands: contlabel       = OFF\n");
  120.          continue;
  121.          }
  122.       if (strcmp(argv[i],"-print") == 0) {        /* printplot */
  123.          printplot  = ON;
  124.          fprintf(stdout,"   Commands: printplot       = ON\n");
  125.          continue;
  126.          }
  127.       if (strcmp(argv[i],"-ps") == 0) {           /* postscript */
  128.          postscript = OFF;
  129.          fprintf(stdout,"   Commands: postscript      = OFF\n");
  130.          continue;
  131.          }
  132.       if (strcmp(argv[i],"-old") == 0) {          /* old format */
  133.          oldformat = ON;
  134.          fprintf(stdout,"   Commands: oldformat       = ON\n");
  135.          continue;
  136.          }
  137.       if (strcmp(argv[i],"-s") == 0) {            /* scale factor */  
  138.          if (++i >= argc) usage(progname);
  139.          scale = atof(argv[i]);
  140.          if (scale < 0.1 || scale > 1) {
  141.             scale = 1.0;
  142.             fprintf(stdout,"   Warning : Invalid scale -");
  143.             fprintf(stdout," scale reset to %f\n",scale);
  144.          }
  145.          fprintf(stdout,"   Commands: scale           = %f\n",scale);
  146.          continue;
  147.          }
  148.       if (strcmp(argv[i],"-o") == 0) {            /* options file */
  149.          if (++i >= argc) usage(progname);
  150.          if ((optionfile =  fopen(argv[i],"r")) == NULL) {
  151.             fprintf(stderr,"cat: can't open %s\n", argv[i]);
  152.             usage(progname);
  153.          } else {
  154.             fprintf(stdout,"   Commands: options-file    = \"%s\"\n",argv[i]);
  155.             parse(optionfile);
  156.             fclose(optionfile);
  157.          }
  158.          continue;
  159.          }
  160.       if (argv[i][0] == '-' && argv[i][1] == 'P') {     /* Printer */
  161.          strcpy(printer,argv[i]); 
  162.          fprintf(stdout,"   Commands: Printer         = %s\n",printer);
  163.          continue;
  164.          }
  165.       if (argv[i][0] == '-') usage(progname);
  166.       if (argv[i][0] != '-' && argv[i][0] != '=') {
  167.          FILEFOUND = TRUE;
  168.          if ((fp = fopen(argv[i], "r")) == NULL) {
  169.             fprintf(stderr,"cat: can't open %s\n", argv[i]);
  170.             exit(1);
  171.          } else {
  172.             imagefile = argv[i];
  173.             fileread(fp);
  174.             fclose(fp);
  175.          }
  176.          }
  177.    }
  178.    if (FILEFOUND == FALSE) usage(progname);
  179.    /* exit(0); */
  180. }
  181.  
  182. /* read in the data */
  183. fileread(fp)
  184. FILE *fp;
  185. {
  186.    float    rdfloat();
  187.    triaptr  insert_tria();
  188.    data     *create_data_array();
  189.    data     get_data_array_value();
  190.    void     free_data_array();
  191.    extern   int     joincurve;
  192.    extern   char    *joinlevel;
  193.    extern   double  xmin, xmax, ymin, ymax, zmin, zmax;
  194.    extern   int     oldformat;
  195.    double   x0,xl,y0,yl,z,delta,ztarget;
  196.    data     *datapt, midpt, pt00, pt01, pt10, pt11;
  197.    int      i,j,nx,ny,isize,jsize; 
  198.  
  199.    zmin = LARGE;
  200.    zmax = -LARGE;
  201.  
  202.    if (!oldformat) {
  203.       /* read the number of x and y points */
  204.       xmin = rdfloat(fp);
  205.       xmax = rdfloat(fp);
  206.       ymin = rdfloat(fp);
  207.       ymax = rdfloat(fp);
  208.       xl   = xmax - xmin;
  209.       yl   = ymax - ymin;
  210.       nx   = (int)rdfloat(fp);
  211.       ny   = (int)rdfloat(fp);
  212.    } else {
  213.       /* old contour format from SPLAT - assume 50 x 50 data */
  214.       x0   = rdfloat(fp);
  215.       y0   = rdfloat(fp);
  216.       xl   = rdfloat(fp);    
  217.       yl   = rdfloat(fp);    
  218.       xmin = x0;
  219.       xmax = x0+xl;
  220.       ymin = y0;
  221.       ymax = y0+yl;
  222.       nx   = 50;
  223.       ny   = 50;
  224.    }
  225.  
  226.    if ((nx < 0) || (ny < 0)) {
  227.       fprintf(stderr,"No of array gridpoints is less than 0 : ");
  228.       fprintf(stderr,"nx=%d ny=%d\n",nx,ny);
  229.       exit(1);
  230.    }
  231.    /* create the new array */
  232.    isize = nx;
  233.    jsize = ny;
  234.    if (joincurve) {
  235.       isize = nx + 2;
  236.       jsize = ny + 2;
  237.    }
  238.    datapt = create_data_array(isize,jsize);
  239.  
  240.    if (joincurve == OFF) {
  241.       /* read the intensity values next */
  242.       /* data goes from i=0 to i=npts-1 */
  243.       for (j=0; j<ny; j++)
  244.       for (i=0; i<nx; i++) {
  245.          z = rdfloat(fp);   
  246.          if (z < zmin) zmin = z;
  247.          if (z > zmax) zmax = z;
  248.          midpt.x = xmin + i*xl/(double)(nx-1);
  249.          midpt.y = ymin + j*yl/(double)(ny-1);
  250.          midpt.z = z;
  251.          assign_data_array_value(datapt,i,j,isize,jsize,&midpt);
  252.       }
  253.       /* now put the values into triangles */
  254.       for (i=0; i<nx-1; i++)
  255.       for (j=0; j<ny-1; j++) {
  256.          pt00 = get_data_array_value(datapt,i  ,j  ,isize,jsize);
  257.          pt10 = get_data_array_value(datapt,i+1,j  ,isize,jsize);
  258.          pt01 = get_data_array_value(datapt,i  ,j+1,isize,jsize);
  259.          pt11 = get_data_array_value(datapt,i+1,j+1,isize,jsize);
  260.          midpt.x = 0.5*(pt00.x + pt10.x);
  261.          midpt.y = 0.5*(pt00.y + pt01.y);
  262.          midpt.z = 0.25*(pt00.z + pt01.z +
  263.                          pt10.z + pt11.z);
  264.          insert_tria(&(pt00),&(pt10),&midpt);
  265.          insert_tria(&(pt10),&(pt11),&midpt);
  266.          insert_tria(&(pt11),&(pt01),&midpt);
  267.          insert_tria(&(pt01),&(pt00),&midpt);
  268.       }
  269.    } else {
  270.       if (xl > yl) delta = 1.0e-4 * yl;
  271.       else         delta = 1.0e-4 * xl;
  272.       /* read the intensity values next */
  273.       /* data goes from i=1 to i=npts   */
  274.       for (j=1; j<=ny; j++)
  275.       for (i=1; i<=nx; i++) {
  276.          z = rdfloat(fp);   
  277.          if (z < zmin) zmin = z;
  278.          if (z > zmax) zmax = z;
  279.          midpt.x = xmin + (i-1)*xl/(double)(nx-1);
  280.          midpt.y = ymin + (j-1)*yl/(double)(ny-1);
  281.          midpt.z = z;
  282.          assign_data_array_value(datapt,i,j,isize,jsize,&midpt);
  283.       }
  284.       if ((strcmp(joinlevel,"HIGH")==0)||(strcmp(joinlevel,"high")==0))
  285.          ztarget = zmax;
  286.       else
  287.          ztarget = zmin;
  288.       /* now fill up a boundary layer, i=0, i=npts+1 */
  289.       for (i=0; i<=nx+1; i++) {
  290.          j=0;
  291.          midpt.x = xmin + (i-1)*xl/(double)(nx-1);
  292.          if (i==0   ) midpt.x = xmin - delta;
  293.          if (i==nx+1) midpt.x = xmax + delta;
  294.          midpt.y = ymin - delta;
  295.          midpt.z = ztarget;
  296.          assign_data_array_value(datapt,i,j,isize,jsize,&midpt);
  297.          j=ny+1;
  298.          midpt.x = xmin + (i-1)*xl/(double)(nx-1);
  299.          if (i==0   ) midpt.x = xmin - delta;
  300.          if (i==nx+1) midpt.x = xmax + delta;
  301.          midpt.y = ymax + delta;
  302.          midpt.z = ztarget;
  303.          assign_data_array_value(datapt,i,j,isize,jsize,&midpt);
  304.       }
  305.       for (j=0; j<=ny+1; j++) {
  306.          i=0;
  307.          midpt.x = xmin - delta;
  308.          midpt.y = ymin + (j-1)*yl/(double)(ny-1);
  309.          if (j==0   ) midpt.y = ymin - delta;
  310.          if (j==ny+1) midpt.y = ymax + delta;
  311.          midpt.z = ztarget;
  312.          assign_data_array_value(datapt,i,j,isize,jsize,&midpt);
  313.          i=nx+1;
  314.          midpt.x = xmax + delta;
  315.          midpt.y = ymin + (j-1)*yl/(double)(ny-1);
  316.          if (j==0   ) midpt.y = ymin - delta;
  317.          if (j==ny+1) midpt.y = ymax + delta;
  318.          midpt.z = ztarget;
  319.          assign_data_array_value(datapt,i,j,isize,jsize,&midpt);
  320.       }
  321.  
  322.       /* now put the values into triangles */
  323.       for (i=0; i<nx+1; i++)
  324.       for (j=0; j<ny+1; j++) {
  325.          pt00 = get_data_array_value(datapt,i  ,j  ,isize,jsize);
  326.          pt10 = get_data_array_value(datapt,i+1,j  ,isize,jsize);
  327.          pt01 = get_data_array_value(datapt,i  ,j+1,isize,jsize);
  328.          pt11 = get_data_array_value(datapt,i+1,j+1,isize,jsize);
  329.          midpt.x = 0.5*(pt00.x + pt10.x);
  330.          midpt.y = 0.5*(pt00.y + pt01.y);
  331.          midpt.z = 0.25*(pt00.z + pt01.z +
  332.                          pt10.z + pt11.z);
  333.          insert_tria(&(pt00),&(pt10),&midpt);
  334.          insert_tria(&(pt10),&(pt11),&midpt);
  335.          insert_tria(&(pt11),&(pt01),&midpt);
  336.          insert_tria(&(pt01),&(pt00),&midpt);
  337.       }
  338.    }
  339.  
  340.    /* now free the array */
  341.    free_data_array(datapt);
  342. }
  343.  
  344. /* Read a float; skip the whole line if failed */
  345. float rdfloat(fp)
  346. FILE *fp;
  347. {
  348.    int n;
  349.    char c;
  350.    float x;
  351.  
  352.    /* try to read a floating variable */
  353.    /* if fail, then skip the whole line */
  354.    while ( (n=fscanf(fp,"%f",&x)) != EOF && n!=1) {
  355.       if (n!=1)
  356.          while ( (c=getc(fp)) != EOF && c != '\n');
  357.       }
  358.  
  359.    if (n==EOF) {
  360.       fprintf(stderr,"Error: Incorrect Format!\n");
  361.       exit(1);
  362.    }
  363.    return(x);
  364. }
  365.  
  366. /* the options */
  367. usage (program)
  368. char *program;
  369. {
  370.       fprintf(stderr,"usage: %s [-o options_file] [-Pprinter] \
  371. [-s scale] [-c level] [-e] [-g] [-j joinlevel] [-l] \
  372. [-print] [-ps] [-old] contour_file\n",program);
  373.       exit(1);
  374. }
  375.